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We present the first determination of the strong coupling constant from the three-jet rate in 
e'^e~ annihilation at LEP, based on a next-to-next-to- leading order (NNLO) perturbative QCD 
prediction. More precisely, we extract aa{Mz) by fitting perturbative QCD predictions at 0{a^s) 
to data from the ALEPH experiment at LEP. Over a large range of the jet-resolution parameter 
j/cut this observable is characterised by small non-perturbative corrections and an excellent stability 
under renormalisation scale variation. We find as{Mz) = 0.1175 ± 0.0020 (exp) ± 0.0015 (theo), 
which is more accurate than the values of as{Mz) from e'^e~ event shape data currently used in 
the world average. 

PACS numbers: 12.38.Bx, 13.66.Bc, 13.66.Jn, 13.87.-a 
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Jet observables in electron-positron annihilation play 
an outstanding role in studying the dynamics of the 
strong interactions described by the theory of Quan- 
tum Chromodynamics (QCD, Q). In particular, jet rates 
and related event-shape observables have been exten- 
sively used for the determination of the QCD coupling 
constant as (see [1, 01 for a review), mostly based on 
data obtained at the e+e" colhders PETRA, LEP and 
SLC at centre-of-mass energies from 14 to 209 GeV. Jets 
are defined using a jet algorithm, which describes how to 
recombine the particles in an event to form the jets. A jet 
algorithm consists of two ingredients: a distance measure 
and a recombination procedure. The distance measure 
is computed for each pair of particles to select the pair 
with the smallest separation in momentum space. If the 
separation is below a pre-defined resolution parameter 
Ucut, the pair are combined according to the recombina- 
tion procedure. The JADE algorithm Q uses the pair 
invariant mass as distance measure. Several improved 
jet algorithms have been proposed for e"'"e~ collisions: 
Durham , Geneva [1] and Cambridge [§] . The Durham 
algorithm has been the most widely used by experiments 
at LEP [10 13] and SLD [TJ], as well as in the reanalysis 
of earlier data at lower energies from JADE [15;]. 

The Durham jet algorithm clusters utilises the distance 
measure 
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for each pair (i,j) of particles, i?vis denotes the energy 
sum of all particles in the final state. The pair with the 
lowest Uij^D is replaced by a pseudo-particle whose four- 
momentum is given by the sum of the four-momenta of 
particles i and j ('i?' recombination scheme). This pro- 
cedure is repeated as long as pairs with invariant mass 
below the predefined resolution parameter Uij^D < J/cut 
are found. Once the clustering is terminated, the re- 



maining (pseudo-)particles are the jets. In experimental 
jet measurements, one studies the jet rates, i.e. jet cross 
sections normalised to the total hadronic cross section, 
as function of the jet- resolution parameter j/cut- 

The theoretical prediction of jet cross sections is made 
within perturbative QCD, where the same jet algorithm 
is applied to the final state partons. The QCD descrip- 
tion of jet production is either based on a fixed-order 
calculation or a parton shower. The fixed order ap- 
proach uses exact parton-level matrix elements includ- 
ing higher order corrections where available and/or an- 
alytical resumniation of large logarithmic corrections for 
a given jet multiplicity. On the other hand, the par- 
ton shower starts with the leading-order matrix element 
for two-jet production and generates higher multiplici- 
ties in an iterative manner, thereby accounting only for 
the leading logarithmic terms from parton-level processes 
with higher multiplicity. In multi-purpose event gen- 
erator programs |16l4l8l |. such parton showers are com- 
plemented by phenomenological models which describe 
the transition from partons to hadrons. These programs 
provide a satisfactory description of multi-jet produc- 
tion rates but, since they generally contain many tunable 
phenomenological parameters, their predictive power is 
limited. Nevertheless, in order to compare parton level 
predictions with experimental hadronic data, these event 
generators are vital to estimate the effects due to hadro- 
nisation and resonance decays. 

Until recently, fixed-order calculations were available 
up to next-to-next-to-leading order (NNLO) for two 
jets Il9l-[2ll and up to next-to- leading order (NLO) for 
three 2a-l24| and four iets [25- 28]. For five and more jets, 
only leading order calculations are available (29I - I3H . For 
jets involving massive quarks, NLO results are available 
for three-jet final states JB^j- The recent calculations of 
the corrections (NNLO) for three-jet production |3a- 
|36| have already led to precise as determinations [37H41| , 
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using event-shape observables measured by ALEPH and 
JADE. However, some of the event-shape variables still 
suffer from a poor convergence of the perturbative ex- 
pansion even at NNLO. Furthermore, the usage of event 
generators, which have been tuned to LEP data, for the 
determination of the hadronisation corrections may lead 
to a bias in the as measurements for some of the event 
shapes [4l[. A comparison of different variables showed 
that jet broadening variables are most affected by missing 
higher orders and a potential hadronisation bias, while 
the differential two-jet rate is most robust against 
these effects, and strongly motivates the present study 
of the three-jet rate. 

In this letter we describe a determination of the strong 
coupling constant from the three-jet rate measured by 
ALEPH nil at LEP. We use the NNLO predictions as 
presented in 3J]. There it was shown that: (i) For large 
values of ?/cut, J/cut > 10~^, the NNLO corrections turn 
out to be very small, while they become substantial for 
medium and low values of t/cut; (ii) The maximum of the 
jet rate is shifted towards higher values of ?/cut compared 
to NLO and is in better agreement with the experimental 
observations; (iii) The theoretical uncertainty is lowered 
considerably compared to NLO, especially in the region 
10~^ > t/cut > 10~^ relevant for precision phenomenol- 
ogy where the theory error is below two per-cent relative 
uncertainty; (iv) Finally, in this t/cut region the parton 
level predictions at NNLO are already very close to the 
experimental measurements, indicating the need for only 
small hadronisation corrections. 

These findings motivate a dedicated analysis of the 
three-jet rate, leading to a precise measurement of as- 
Our analysis closely follows the procedure described in 
[H El. The ALEPH data [13] at LEP are based on 
the reconstructed momenta and energies of charged and 
neutral particles. The measurements have been corrected 
for detector effects, i.e. the final distributions correspond 
to the so-called particle (or hadron) level, and for initial 
state photonic radiation. In the simulation of the de- 
tector response to particles, a bias is introduced by the 
choice of the physics event generator. This leads to a sys- 
tematic uncertainty on the three-jet rate of about 1.5% 
for the relevant t/cut range. Further experimental sys- 
tematic effects are estimated by a variation of the track- 
and event-selection cuts as advocated in [42:], giving an 
additional small systematic uncertainty of about 1%. 

We construct the perturbative expansion up to ©(a^) 
as described in [4l|, with the coefficients obtained from 
[S^ l- These are valid for massless quarks. We take into 
account bottom mass effects up to NLO ;32|], for a pole 
b-quark mass of Afb = 4.5 GeV. The latter is varied by 
±0.5 GeV in order to estimate the impact of the b-quark 
mass uncertainty on the value of the strong coupling. For 
the normalisation to the total hadronic cross section CThad 
we follow the procedure adopted in [Zlll, which is based 
on a N^LO calculation (O(a^) in QCD) for tThad [H, 
including mass corrections for the b-quark up to 0{as) 
and the leading mass terms to 0{al). Weak corrections 



to the three-jet rate were computed very recently [44|. 
They are at the one per-mille level for Q = Mz and are 
neglected here. 

The nominal value for the renormalisation scale Xf^ — 
/i/Q is unity. It is varied between 0.5 < a;^ < 2 in order 
to assess the systematic uncertainty related to yet un- 
known higher order corrections. No attempt is made to 
combine the NNLO predictions with resummation calcu- 
lations. At present, the resummation of the three-jet 
rate )Ajs only fully consistent at leading logarithmic 
level (45|, and resummation effects only become numeri- 
cally relevant over fixed-order NNLO for In t/cut ^ — 4.5 
(as can be seen from the Ya transition parameter distri- 
bution [1^), which is below our region of interest. 

In order to compare the perturbative parton level 
thoretical prediction with the hadronic data, it is nec- 
essary to apply a correction for hadronisation and res- 
onance decays. This bin-by-bin correction is computed 
with the PYTHIA JJ], HERWIG [17] and ARIADNE 
(isj Monte Carlo generators, all tuned to global hadronic 
observables at Mz [13] • The parton level is defined by 
the quarks and gluons present at the end of the parton 
shower in PYTHIA and HERWIG and the partons re- 
sulting from the colour dipole radiation in ARIADNE. 
Our central values for the strong coupling constant are 
obtained with hadronisation corrections from PYTHIA, 
which are at the level of 5%. We define the systematic 
uncertainty on as{Mz) due to these hadronisation cor- 
rections as the biggest deviation observed when using any 
of the other generators. Motivated by the observations 
in [ill, we verified that the shapes of the Monte Carlo 
parton level predictions are in fair agreement with those 
at NNLO, for reasonable choices of the strong coupling. 
Furthermore, the ratios of these predictions are relatively 
flat over the relevant t/cut range, giving further confldence 
in the reliability of the hadronisation corrections. 

The corrected ALEPH measurements for the three-jet 
rate are compared to the theoretical calculation at par- 
ticle level. Values for as{Mz) are obtained by a least- 
squares fit, performed separately for each t/cut value in the 
range listed in Table U (for the data at the Z peak), to- 
gether with the uncertainties as described above. These 
results are also displayed in Fig. [T] We observe a nice 
stability of the results, within their total uncertainties, 
down to resolution parameters of In t/cut ~ —4.5. Beyond 
that value we find a fall-off of as{Mz), most likely re- 
lated to the onset of large logarithmic corrections from 
higher perturbative orders, which are not accounted for 
in our perturbative prediction. 

As final result we quote our measurement for t/cut — 
0.02, which represents an optimal compromise between 
minimal systematic uncertainty and stability. We find 

as{Mz) = 0.1175 ± 0.0020 (exp) ± 0.0015 (theo) 

where the first uncertainty includes (in quadrature) the 
contributions from statistics, detector corrections and 
experimental selection cuts, and the second error is 
the quadratic sum of b-quark mass and renormalisa- 
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FIG. 1: Determinations of as{Mz) from the three-jet rate, 
measured by ALEPH at the Z peak, for several values of 
the jet-resolution parameter j/cut- The error bars show the 
statistical uncertainty, whereas the shaded band indicates the 
total error, including the systematic uncertainty. The various 
contributions to the latter are displayed in the lower plot. 



tion scale uncertainties (of. Table |T|. We also per- 
formed similar measurements for the LEP2 energies be- 
tween 133 and 206 GeV, where we find consistent val- 
ues for as{Mz), but with considerably larger statisti- 
cal uncertainties. Combining the errors in quadrature, 
yields as{Mz) = 0.1175 ± 0.0025 which is in excellent 
agreement with the latest world average value [41 of 
as{Mz) = 0.1184 ± 0.0007 that is based on a number 
of measurements from r-decay, lattice gauge theory, Up- 
silon decay, DIS and e+e" data. As expected, our the- 
oretical uncertainty is smaller than that obtained from 
fits of event-shape distributions, and even smaller than 
the experimental error, which is dominated by the model- 
dependence of the detector corrections. Our result is also 
more precise than the two extractions of from e^e~ 
event-shape data [13, |4l| currently used in the world av- 
erage 

In this letter we reported on the first determination 
of the strong coupling constant from the three-jet rate 



in e~^e~ annihilation at LEP, based on a NNLO per- 
turbative QCD prediction. We find a precise value of 
as(Mz) with an uncertainty of 2%, consistent with the 
world average. This verifies the expectations that the 
three-jet rate is an excellent observable for this kind of 
analysis, thanks to the good behaviour of its perturbative 
and non-perturbative contributions over a sizable range 
of jet-resolution parameters. 
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-5.1 
-4.9 
-4.7 
-4.5 
-4.3 
-4.1 
-3.9 
-3.7 
-3.5 
-3.3 
-3.1 
-2.9 
-2.7 
-2.5 
-2.3 
-2.1 
-1.9 
-1.7 
-1.5 
-1.3 



:{Mz) stat. dot. cxp. 



had. 



mass 



pert. total 



0.1110 
0.1124 
0.1147 
0.1153 
0.1159 
0.1170 
0.1175 
0.1179 
0.1183 
0.1184 
0.1179 
0.1177 
0.1180 
0.1169 
0.1166 
0.1166 
0.1191 
0.1173 
0.1175 
0.1159 



0.0004 
0.0004 
0.0004 
0.0004 
0.0004 
0.0004 
0.0004 
0.0004 
0.0004 
0.0004 
0.0004 
0.0004 
0.0004 
0.0005 
0.0005 
0.0006 
0.0008 
0.0010 
0.0016 
0.0037 



0.0013 
0.0015 
0.0015 
0.0015 
0.0016 
0.0016 
0.0016 
0.0016 
0.0015 
0.0015 
0.0016 
0.0019 
0.0020 
0.0021 
0.0019 
0.0020 
0.0021 
0.0015 
0.0005 
0.0014 



0.0008 
0.0007 
0.0008 
0.0008 
0.0009 
0.0009 
0.0011 
0.0011 
0.0009 
0.0011 
0.0013 
0.0013 
0.0013 
0.0015 
0.0018 
0.0020 
0.0019 
0.0023 
0.0029 
0.0029 



0.0003 
0.0003 
0.0004 
0.0005 
0.0005 
0.0005 
0.0006 
0.0006 
0.0006 
0.0008 
0.0010 
0.0010 
0.0013 
0.0013 
0.0014 
0.0015 
0.0014 
0.0016 
0.0014 
0.0018 



0.0004 
0.0003 
0.0003 
0.0003 
0.0003 
0.0003 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0001 
0.0001 
0.0001 
0.0001 
0.0002 
0.0001 
0.0001 
0.0004 



0.0020 
0.0013 
0.0012 
0.0006 
0.0010 
0.0012 
0.0014 
0.0016 
0.0018 
0.0019 
0.0021 
0.0021 
0.0020 
0.0021 
0.0021 
0.0020 
0.0016 
0.0019 
0.0017 
0.0011 



0.0025 
0.0022 
0.0022 
0.0019 
0.0022 
0.0023 
0.0025 
0.0026 
0.0026 
0.0029 
0.0031 
0.0033 
0.0034 
0.0036 
0.0037 
0.0038 
0.0036 
0.0038 
0.0040 
0.0054 



TABLE I: Results of Qs {Mz) extracted from the three-jet rate 
measured by ALEPH at LEPl. The uncertainty contribu- 
tions are given for the statistical error (stat.), the uncertainty 
related to the choice of the generator for the simulation of 
the detector response (det.), the quadratic sum of all other 
experimental systematic uncertainties arising from track and 
event selection cut variations (exp.), the hadronisation un- 
certainty obtained by the maximum difference between either 
PYTHIA, HERWIG or ARIADNE (had.), the uncertainty on 
the b-quark mass correction procedure (mass) and the un- 
certainty for missing higher orders (pert.) estimated by a 
variation of the renormalisation scale. 
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